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Abstract 



Error performance is one of the main performance measures and derivation of its closed-form 
expression has proved to be quite involved for certain systems. In this letter, a unified closed-form 
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expression, applicable to different binary modulation schemes, for the bit error rate of dual-branch 
selection diversity based systems undergoing independent but not necessarily identically distributed 
generalized-^ fading is derived in terms of the extended generalized bivariate Meijer G-function. 

Index Terms 

Diversity schemes, selection combining, dual-branch selection diversity, binary modulation schemes, 
generalized-^T (GK) model, composite fading, bit error rate (BER), and Meijer G-function distribution. 

I. Introduction 

In recent times, different diversity schemes have taken up an important position in the wireless 
communication systems. The main reason behind this is that these different diversity schemes 
allow for multiple transmission and/or reception paths for the same signal [1]. One of the simplest 
diversity combining scheme is the selection combining (SC) diversity scheme where only one 
of the diversity branches is processed. Specifically, SC scheme chooses the branch with highest 
signal-to-noise ratio (SNR) [2]-[3]. 

Additionally, wireless communications are driven by a complicated phenomenon known as 
radio-wave propagation that is characterized by various effects including multipath fading and 
shadowing. The statistical behavior of these effects is described by different models depending 
on the nature of the communication environment. It is becoming necessary to study such 
effects i.e. large-scale fading as well as small-scale fading concurrently as the multihop relay 
networks are emerging in the current times. Using the Nakagami multipath fading model that 
is versatile enough to model various multipath fading conditions ranging from severe fading to 
non-fading scenario, and the Gamma model for shadowing has led to the generalized-^ (Gamma- 
Gamma) composite fading model [4]-[8]. Generalized-^ (GK) distribution, earlier used in radar 
applications and recently being used in the context of wireless digital communications over fading 
channels, is one of the relatively new tractable models used to describe the statistical behavior 
of multipath fading and shadowing effects as compared to log-normal based models. GK fading 
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model is quite general model as it includes ^T-distribution as its special case and accurately 
approximates many other fading models such as Nakagami-m and Rayleigh-Lognormal (R-L) 
([6] and references therein). Finally, GK distribution is a distribution of the product of two 
independent Gamma random variables (RV) and hence is a special case of the Fox //-function 
and in turn a special case of Meijer G-function where the product of two Meijer G-function's 
can be represented in terms of a extended generalized bivariate Meijer G-function (EGBMGF) 
(refer Table I) [9]. 

It is noteworthy to mention that bit error rate (BER) is one of the most important performance 
measures that forms the basis in designing wireless communication systems. Based on the open 
technical literature and upto the best of our knowledge, error analysis has been performed for 
dual diversity with SC over log-normal fading channels in closed-form using moment generating 
function (MGF) based approach in [10] and with Weibull fading channel as an approximate using 
characteristic function (CF) based approach in [11]. Additionally, error performance analysis 
of SC systems with independent and identically distributed (i.i.d.) GK fading branches was 
performed in [12] involving integral form expressions. Further, in [13] the analysis was performed 
for dual-branch SC citing the difficulty in deriving the expression for the probability density 
function (PDF). This issue was tackled in [14] for an arbitrary number of branches and the 
authors therein have described and utilized a method to perform the BER analysis directly from 
the cumulative density function (CDF) eliminating the need of deriving the PDF and relying on 
the Gauss-Laguerre quadrature technique. 

In this work, we revisit this problem under the umbrella of the iZ-functions and derive exact 
closed-form expression of the BER of binary modulation systems with dual-branch SC scheme 
and undergoing GK fading where the channels are independent but not necessarily identically 
distributed (i.n.i.d.). The remainder of the paper is organized as follows. Section II introduces 
the system and the GK channel model. Next, section III presents some statistical characteristics 
of GK fading channel model followed by the analytical BER analysis in section IV, and finally, 
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section V discusses the results followed by the summary of the paper in section VI. 

II. The Generalized-^ Fading System and Channel Model 

A SC based communication system with a source and a destination is considered with i.n.i.d. 
channels as follows 

Y = aX + n, (1) 

where Y is the received signal at the receiver end, X is the transmitted signal, a is the channel 
gain, and n is the additive white Gaussian noise (AWGN). In a Nakagami multipath fading 
channel, 7 = \a\ 2 follows Gamma distribution; additionally, the shadowing component is also 
assumed to follow a Gamma distribution. Hence, the channel gains experience composite fading 
whose statistics follow a generalized-A' distribution given by 



2l) m m+rn, 



PM = rvtl ^ , 7 an?!Ba " 1 ^m.- mm (26V7). ( 2 ) 



T{m m )T{m s 



where T(-) is the Gamma function as defined in [15, Eq. (8.310)], m m and m s are the Nakagami 
multipath fading and shadowing parameters, respectively. In (2) K m (-) is the modified Bessel 



function of the second kind and order m,b— y m ^" s , and f2 is the mean of the local power. The 
parameters m m and m s quantify the severity of multipath fading and shadowing, respectively, 
in the sense that small values of m m and m s indicate severe multipath fading and shadowing 
conditions respectively, and vice versa. The instantaneous SNR of the n th branch is given by 
7„ = (Eb/No) %n where x n is the signal amplitude for the n th branch, Ef, is the average energy 
per bit and N is the power spectral density of the AWGN. 

III. Statistical Characteristics 

The PDF and CDF expressions of the GK RVs can also be written in terms of Meijer G- 
function. 

Lemma 1: The PDF of a GK RV can be expressed in terms of Meijer-C7 function as 

'm m m s 



m m m s \ 2,0 



fin 



V 



m m -l,m s -\ 



y>0. (3) 
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Proof: We may use the fact that the PDF of the product of iV independent Gamma RVs 



can expressed as a //-function PDF that is given by [16, Eq. (6.4.9)] as follows 



n i 



X 



(fci-1,1), ...,(^-1,1) 



x > 0. 



(4) 



Then, with N = 2, the GK PDF can be expressed as 



p(x) 



Tn m Tn s 



T(m m )T(m s )n J °' 2 



H, 



2.0 



TYl m TYl s 
^0 



(m m - 1, 1), (m s - 1, 1) 



x > 0. (5) 



Now by applying [16, Eq. (6.2.8)], the expression in (3) follows. ■ 
Further, substituting (3) in [17, Eq. (26)] and utilizing [16, Eq. (6.2.4)], the CDF of GK can 
be written as 

1 



1 



r (m m ) r (m a 



^1,3 



m m m s 
fin 



7 



m m ,m s ,0 



,7>0, 



(6) 



where G[-] is the Meijer G-function [15]. 



IV. BER Analysis 

In SC combining scheme, the highest SNR branch is selected. In our case, for dual-diversity, 
the SNR, 7 SC , is given by 

7 SC = max(7i,7 2 ). (7) 

The CDF of 7 SC is given by 

2 

F( 7sc ) = Pr(max( 7l ,7 2 ) < 7sc ) = Y[F^(j ac ). (8) 
The BER for SC is given by 



n=l 



Pe (e|7«.) / 7 „ (7) d-f sc = Pe (e|7sc) dF ln {^ sc ) 
^0 



(9) 



where P e (e|7 ac ) is the conditional error probability (CEP) for the given SNR. A unified CEP 
expression for coherent and non-coherent binary modulation schemes over an AWGN channel 
is given in [18] as 

P r £ U) = JWi), (10) 



2F(p) 
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where ■) is the complementary incomplete gamma function [15, Eq. (8.350.2)]. The parame- 
ters p and q in (10) account for different modulation schemes. For an extensive list of modulation 
schemes represented by these parameters, one may look into [19]. Now, applying integration by 
parts in (9), we get 



POO 

Pe = Pe (e| 7s c) F(lsc)\™ ~ / F(lsc)dP e (e| 7ac ) . 

J 



(11) 



The first term goes to zero using [20, Eq. (6.5.3)]. Further, substituting (10) into (11) and using 
[20, Eq. (6.5.25)], the average BER can be written as 



exp(-g7 sc )7^ 1 F(j sc )d% 



2r(p) Jo 

On substituting (8) in the above obtained expression, we get 



(12) 



P P 



q 



2T(p) J 



exp(-g 7sc ) 7 f c 1 Yl F ln {pf) d lsc- 



(13) 



n=l 



Using [9], we obtain the product of the CDFs present in the above expression in terms of 
EGBMGF as 



\F ln {l) = ^(7)^(7) = «i S 



n=l 





0,0 
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(14) 



where S[] is the EGBMGF as given in [21, Eq. (2.1)], k x 



1 



m mi ,m si ,0, k 3 = m m2 ,m S2 ,0, k 4 = ""T^ 1 and k 5 
be expressed as, 



r(m mi )r(m si )r(m m2 )r(m S2 ) ' 
"^ ms2 • The above expression can also 
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where G(-) is EGBMGF as in [22]. Additionally, (14) or (15) can be represented as, 



n=l 



0,0 
0,0 



'2,1^ 

V 1 ' 3 / 



1 

K 2 



where S[] is EGBMGF as in [23, Eq. (4)]. 



V 1 ' 3 / 



^3 

(16) 



Lemma 2 [21, Eq. (2.1)]: The integral involving the EGBMGF of two variables with an 
exponential term with the RV as one of its argument and a term with RV itself evaluates to 

p,0 
A-p,B 



\-i e -»x s 



(27T) 



J(i-p)£ 
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(17) 



where Afa A) = ± ■ 

Now, substituting (14) or (15) or (16) into (13), then using the Lemma 2 given above and 
performing additional manipulations, we get the desired closed-form expression for the average 
BER as 
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(18) 
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or equivalently 
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1 
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V. Results and Discussion 

The numerical results for BER of SC scheme with dual-diversity over i.n.i.d. GK fading 
channels are presented in this section. 

The exact solution presented above in (18), (19) and/or (20) has not been found computable and 
hence its computability/evaluation was implemented using Mathematica as can be seen in Table 
II. With this implementation, the EGBMGF can be evaluated fast and accurately. An illustrative 
code is shown in Table III. This computability, therefore, has been utilized for different digital 
modulation schemes and is employed to discuss the results in comparison to respective Monte 
Carlo simulation outcomes. 

The average SNR per bit in all the scenarios discussed is assumed to be equal. In addition, 
different digital modulation schemes are represented based on the values of p and q where 
p = 0.5 and q — 1 represents binary phase shift keying (BPSK), p — 1 and q — 1 represents 
differential phase shift keying (DPSK) and binary frequency shift keying (BFSK) is represented 
by p = 0.5 and q = 0.5. In Monte Carlo simulations, the GK fading channel was generated by 
the product of two independent gamma RVs. 

We observe from figure 1 that this implemented computability of EGBMGF provides a perfect 
match to the MATLAB simulated results and the results are as expected i.e. the BER increases 
as the shadowing effect increases (i.e. value of m s decreases) while keeping multipath fading 
constant at m m = 1. The figure shown represents BPSK. Its important to note here that these 
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values for the parameters were selected randomly to prove the validity of the obtained results 
and hence specific values based on the standards can be used to obtain the required results. 

Similar outcomes may be obtained for BFSK and DPSK. Additionally, similar analysis can 
be done for constant shadowing effect and varying the multipath fading. 

Furthermore, to demonstrate the case that the results presented here also handle the presence 
of i.n.i.d. GK channels, following figure 2 presents the different modulation schemes with 
different effects of multipath fading and shadowing on both their channels. The values utilized 
for multipath fading and shadowing were as follows; m mi = 1, m m2 = 2, m si = 0.5, and 
m S2 = 4. It can be seen that, as expected, BPSK outperforms the other modulation schemes and 
BFSK and DPSK perform in similar fashion at lower SNR whereas as the SNR increases DPSK 
performs better than BFSK. 

Similar results for any other values of m' m s and m' s s can be observed for the exact closed-form 
BER for dual-diversity i.n.i.d. GK channels presented in this work. 

VI. Concluding Remarks 

An exact closed-form expression for the BER performance of different binary modulations 
with dual-branch SC scheme over i.n.i.d. GK fading was derived. The analytical calculations 
were done utilizing a general class of special functions, specifically, the EGBMGF. In addition, 
this work presents numerical examples to illustrate the mathematical formulation developed in 
this work and to show the effect of the fading and shadowing severity and unbalance on the 
system performance. 
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TABLE I 

Extended Generalized Bivariate Meijer G-Function (EGBMGF) 



Representation 1: Based on [21] 




(a); (6) 
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where A + C"<S + D, ,4 + £<B + F. 



Representation 2: Based on [23] 











mi, 






( 








( 






X 


= s 






a pi 






C P2 




n 3 , m 3 1 


r " 






















y 








Pi, 9i 


b qi 




\ 


f>2,?2 y 


d q2 




\ 


f>3,<73 y 


J 93 



G 



mi, 0:712, ^2:^3,^3 
Pl,9i:P2,92:P3>93 



«1. 

6i, 



Cl,. 

di, 



ei, 
/i, 



■ , °P3 

■ , /93 
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Wit r tfj - t)x'y*dsdt 



n^ ro3+ i r fe-*)n*L„ 3+ ir(i -h+ty 



where C\ and C2 are two suitable contours and positive integers pi, P2, P3, ?i, <?2, <?3, itoi, m?, m?,, 712, and 713 satisfy the 
following inequalities. <?2 > 1, 93 > 1, pi > 0, < mi < pi, < rri2 < P2, < 712 < §2, < 7713 < P3, < 713 < 93, 
Pi + P2 < <?i + 92, Pi + P3 < <Zi + 93- The values a; = and y — are excluded. 

It may be useful to learn the relationship between both the representations shown above. Following equalities must be noted 
from both the above representations; p = mi, A — pi, B — q\, q — ma, r = 712, C = P2, D = q2, k = 7713, Z = 713, i? = P3, 
and F = qz. 
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TABLE II 

Extended Generalized Bivariate Meijer G-Function (EGBMGF) Mathematica Implementation 



(♦Extended Generalized Bivariate Meijer G-Function (EGBMGF) *) 
Clear All; 
(♦Exception*) 

S : : InconsistentCoef f s = "Inconsistent coefficients!"; 

S[{ast_, bst_} , {as_, bs_} , {at_, bt_} , {zs_, zt_}] :=Module[{}, 

(* Gamma product terms with only 's' as argument with other parameters *) 
Pas = Function [u, Product [Gamma [1 - as [ [1 , n]] + u], {n, 1, Length [as [ [1] ] ] } ] ] ; 
Qas = Function [u, Product [Gamma [as [ [2 , n] ] - u], {n, 1 , Length [as [ [2] ] ] } ] ] ; 
Pbs = Function[u, Product [Gamma [bs [ [1, n] ] - u ] , {n, 1, Length[bs[ [1] ] ] }] ] ; 
Qbs = Function [u, Product [Gamma [1 -bs[[2, n]] + u], {n, 1, Length [bs [ [2] ] ] }] ] ; 
Ms = Function [u, Pas [u] Pbs[u] / (Qas [u] Qbs[u])] ; 

(*Gamma product terms with only 't' as argument with other parameters *) 
Pat = Function [u, Product [Gamma [1 - at[ [1, n] ] + u ] , {n, 1, Length [at [[ 1] ]]}]]; 
Qat = Function [u, Product [Gamma [at [ [2 , n] ] - u], {n, 1, Length [at [ [2] ] ] } ] ] ; 
Pbt = Function[u, Product [Gamma [bt [[ 1 , n] ] - u ] , {n, 1, Length [bt [[ 1] ]]}]]; 
Qbt = Function[u, Product [Gamma [1 - bt[ [2, n] ] + u ] , {n, 1, Length [bt [ [2] ]]}]]; 
Mt = Function [u, Pat[u] Pbt[u] / (Qat[u] Qbt[u])] ; 

(*Gamma product terms with 's+t' as argument with other parameters *) 
Past = Function [u, Product [Gamma [ast [[ 1 , n] ]+ u] , {n, 1 , Length [ast [ [1] ]]}]] ; 
Qast = Function [u, Product [Gamma [ 1 - ast [ [2 , n] ] - u], {n, 1, Length [ast [ [2] ]]}]] ; 
Qbst = Function [u, Product [Gamma [bst [ [2 , n] ] + u], {n, 1, Length [bst [ [2] ]]}]] ; 
Mst = Function [u, Past[u] / (Qast[u] Qbst[u])] ; 

(♦Countour limiters (Depends on numerator Gamma 

arguments i.e. it must be half of the least valued Gamma arguments)*) 
Rs = 1 / 4 ; 
Rt = 1 / 4 ; 

( ^Assignments and Declarations*) 
Zs = zs; 
Zt = zt; 
W = 50; 

(*Final Evaluation*) 
Print [ "Numerical Integration : " ] ; 
1 

value = NIntegrate[MT[s, t] Zs s Zt b , {s, Rs - I W, Rs + IW}, {t, Rt-IW, Rt + I W}] ; 

(2ttI) 2 

(♦Returning back the value*) 
Return [value] ; 

]'• 

(*End of EGBMGF*) 
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TABLE III 

Extended Generalized Bivariate Meijer G-Function (EGBMGF) Mathematica Implementation Usage 

Example 



Prior moving to the illustration, it may be useful to know the relationship between the representations shown in Table I and the 
implemented code presented in Table II. Following equalities must be noted from both the representations and the implemented 
code; (a) = a Pl = ast_, (b) — b qi = bst_, (c) = c P2 = as_, (d) = d q2 = bs_, (e) = e P3 = at_, and (/) = f Q3 = bt_. 



(♦Testing*) 
( *Declarations*) 
p = . 5 ; q = 1 ; 

mml = 1; msl = 2; mm2 = 1; ms2 = 2 ; 
Dl = 1; £22 = 1; 

snr = 10 A (15 / 10) (*SNR =0-20 dBs*) ; 
(♦Invoking the implemented EGBMGF module*) 



An illustration of the EGBMGF implemented code usage 



1 



s[{{{p}, {}}, {{}, {}}}, 



B = 



2 Gamma [p] Gamma [mml] Gamma [msl] Gamma[mm2] Gamma[ms2] 
{{{1}, {}}, {{mml, msl}, {0}}}, 
{{{1}, {}}, {{mm2, ms2}, {0}}}, 




mml msl 



mm2 ms2 



(*END*) 



0.00102393 - 7.09829x10 
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Fig. 1. I.I.D. BPSK BER for m m = 1 and varying m s . 
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Fig. 2. BER for different modulation schemes undergoing i.n.i.d. channels with m mi = 1, m m2 — 2, m si = 0.5, and 
m S2 = 4. 
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